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Abstract 

In this paper, we present the point symmetry group of three-dimensional 
homogeneous Helmholtz equation, when we consider the cylindrical coor- 
dinate system. In continuation, we present a complete set of functionally 
independent invariants of the equation along with the form of the general 
solution provided by these invariants. Finally, we find an optimal system 
of one-dimensional Lie subgroups of the full symmetry group. 
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1 Introduction 

The Helmholtz equation, named for Hermann von Helmholtz, is the elliptic 
partial differential equation (in homogenous case) 

Vit + fc^u 0, 

when u is a vector function and V is the vector Laplacian. 

The three-dimensional homogeneous cylindrical Helmholtz equation (in brief 
CHE), is a Helmholtz differential equation defined in the cylindrical {r,6,z) 
coordinate system, i.e. 

CHE : Urr + -Ur + —Ugg + U^z + = 0, (1) 

where u is a function of variables r, 9, z and fc is a constant. 

The Helmholtz equation naturally appears from general conservation laws 
of physics and often arises in the study of physical problems involving partial 
differential equations (PDEs) in both space and time, more particularly in me- 
chanics, acoustics, electromagnetics, etc. (see for details 0111111]). 
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This equation can be interpreted as a wave equation for monochromatic 
waves and can be derived from the heat conduction equation, Schrodinger equa- 
tion, telegraph and other wave type, or evolutionary equations [3]. When k — 0, 
the Helmholtz differential equation reduces to Laplace's equation and when 
fc^ < (for imaginary), the equation becomes the space part of the diffusion 
equation. 

In the mathematical sense, it appears as an eigenvalue problem for the 
Laplace operator (see [31 IH H]). As is well-known, there are some math- 
ematical theory which can be applied to a solution of the Helmholtz equation in 
Cartesian, spherical, cylindrical coordinates, and etc., by applying the technique 
of separation of variables and boundary element method, its various solutions 
for many basic shapes, for example in a vibrating membrane, along with other 
different methods [l] [H |4l [9] . In addition, since the Helmholtz equation is an 
elliptic equation, so it is usual to consider boundary value problem along with 
transmission conditions and integral representation of solutions (see O Sj and 
references therein). 

At the present work, we are dealing with the symmetry analysis and a one- 
dimensional optimal system of the Helmholtz equation in cylindrical coordinate. 
Cylindrical coordinates are the natural choice for studying different shapes of 
the Helmholtz equation, that is a good motivation for work with this coordinate. 
Also we find a form of general solutions providing by the differential invariants 
of the equation. 



In this section, we are deal with symmetry analysis of CHE, when we act the 
point symmetry group on the second order jet space = J-^(R'^,]R) with local 
coordinates (r, 9, z, u, u'-^', it^^)). In these charts m*^"^ (a = 1, 2) indicates the set 
of all a*^ derivatives of the dependent variable u with respect to independent 
variable r,9, z. 

The base space of the jet bundle is the space Ai = X xU independent 
and dependent variables (a;, u) for x = (r, 0, z), that is in fact the zero-order jet 
space J°(R'^, K). CHE defines a subvariety of jet bundle of three-dimensional 
submanifold of M, that is, graph of function u — u{r, 9, z). 

The point symmetry group on Ai is defined by transformations in the general 
form of 

r — ip{r,9, z,u), 9 ^ tlj{r,9, z,u), z — xir,9, z,u), u — ({r,9, z,u), 
for arbitrary smooth functions ip^ip^x ^-^d C- 



be the general form of an infinitesimal generator of point symmetry transforma- 
tions on A^, in which the coefficients are arbitrary smooth functions. Let the 



2 Symmetry of CHE 
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second order prolongation of v be defined as 

where we assumed that each index signifies to the derivation in respect to the 
index, and J = (ji, ^2) is a midti-index with length < #J < 2. In addition, 
coefficients are given by the prolongation formula (see [lOj ) 



i=l 



In the latter relation, Dj is the total derivative with respect to the variables 
signified by J, and Q = i] — X]i=i is the characteristic of v. 

According to [lOj , the vector field v is an infinitesimal symmetry of Eq. ([1]) , 
if and only if it satisfies the infinitesimal invariance condition 



^Urr + -Ur + —Ugg + U^z + k^U^ — 0, (3) 



whenever Eq. ^ vanishes. 

Eq. ([3]) leads to the infinitesimal determining equation, that must be solved 
for arbitrary functions and ry, when they are dependent to variables r, 0, z, u 
and not to the derivatives u in respect to its variables. Hence, the induced 
equation satisfies for and 77, if and only if the following determining equations 
are fulfilled 

£^ =0 £^ =0 f ^ = £^ =0 

a^O, Cf. = 0, £1 = 0, 

r'£i+£'e = o, ei = o, ^1 = 0, e.+er = o, 

2 + rry„„ - 2 r^^^ = 0, r^. - r£l - = 0, " ^1 - 0, 
2eL-?/u« = 0, r^l+£l=Q, ?™ + eL = 0, (4) 
Ceu + ^"^^tu = Oj ~ 2 ^gjj = 0, ''^Cru + Ce« = 0' 

r2 k^u £l + 2rjeu- - ^le " r'C " ^'CL = 0, 

3 r2 fc2^, ^1 + ^ ^1 + 2 r^ry™ - r^e.z " " Cee " '^'C.V = 0, 

k'u£l + 2r2ry,„ - r£'^ - r^Q. ~ ^ee - r'C = «, 
r^Vz.z + 2 " ^u) + ''^^r-r + rrjr + 1100 + r^ri k"^ = 0. 

The solution space g of the system of determining equations (jl]) with respect 
to the unknown coefficients of v, is the Lie algebra of infinitesimal symmetry of 
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Eq. ([T|). By solving system (jH) (that is applicable with Maple) we obtain the 
below solutions 

= (ci z + C3) sin0 + (c2 z + C4) COS0, 

= -i(c2 z + C4) sin0 + i(ci z + C3) cos6' + C5, 
r r 

= —r{ci sinO + C2 cos9) + cq, (5) 
77 = exp(-^ci'z - Vc2^')(c7exp(2^/ci'e') + C8)(c9exp(2i/ci'z) + cio) 

X (cii BY(V-C2, + fc^r) + C12 BJ(V-C2, ^03 + fc^r)^ + C13 u, 

where, BJ and BY are the Bessel functions of the first and second kinds, resp. 
which satisfy Bessel's equation x'^ y" + xy' + (x^ — v'^)y = 0, and Ci (for 1 < i < 
13) are arbitrary constants. 

Noting to the concluded form of infinitesimal operators of Eq. ([T]), one may 
divide the resulting form of v to the set of following vector fields, as a basis of 
Lie algebra g of the symmetry group G of CHE 

Y ^ Y ^ Y ^ 

0& oz ou 

. ^ d cos 6 d ^ d sm9 d 

X4 = sm fc* — H — , X5 = - cos fc* — H — , 

or r of) or r aO 

^ d z smO d ^ d 

Xe = -z cose — H — + r cos6' — , 

or r OO oz 

. „ d z cos 6 d ■ ^ d 

At ^ z smO — -\ — ~ r sm& —. 

Or r oO oz 

The commutator of every Xi and X.j i.e. [Xi, Xj] :— Xi Xj — Xj Xi (1 < i,j < 
7), is straightforwardly a linear combination of vector fields (O. The complete 
set of commutators is given in Table [T] 





Xi 


X2 


X3 


Xi 


X5 


Xe 


X7 


Xi 











-X5 


Xi 


X7 


-Xe 


X2 

















X5 


Xi 



























^5 

















-X2 


X5 


-Xi 














-X2 





Xe 


-Xj 


-X5 








X2 





Xi 


X7 


Xe 


-Xi 





X2 





-Xi 






Table 1: Commutator table for g. 
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Theorem 1. The complete set of infinitesimal generators of point symmetry 
group of Eq. {Ip, is introduced in (0) and ^ and is a seven- dimensional Lie 
algebra with commutators given in Table\^ 

Since the derived algebra g^^^ = [g, g] = q/{X^) and for each k > 1, gC^'+i) = 
= so the descending sequence of the derived subalgebras of g, 

030(1)^0(2)^... 

does not terminate with a nuU ideal and hence is not solvable [51 [T^] . Also, 
is not semisimple, since by applying Cartan's criterion of semisimplicity [12], 
its Killing form defined for v = '^^ViXi and w — 'Y^- WjXj in as 

K{v, w) — tr(ad(u) o ad(?i;)) — — 4(uiU'i + vqWq + viwi) 

is degenerate. 

The radical of is equal to the Lie subalgebra r := (X2, X3, X4, X5) and 
hence Levi decomposition of is in the form of = r s, for semisimple 
subalgebra 5 := q/z — (Xi, Xq, Xj). According to the Table [1] r is an abelian 
group of dimension 4, thus is isomorphic to B.'^, and the (real) Lie subalgebras is 
isomorphic to the Lie algebra so (3) of the special orthogonal group. In addition, 
s is not an ideal of and therefore we find the structure of the symmetry algebra 
of CHE as following semi-direct product 







X so(3). 



(7) 



For finding the symmetry group G of the Lie algebra generated by infinitesimal 
operators (O, we determine the one-parameter group (flow) of these vector 
fields: 



Theorem 2. If we denote the group transformations of each infinitesimal gen- 
erator Xi (i — 1, - ■ ■ ,7) with gi{s) for different values of parameter s, then we 
have 



92{s) 
93{s) 



(r, 6, z, u) I — > (r, 9 s, z,u), 
(r, 9, z, u) I — > (r, 9,z + s, u), 
(r, 9, z, u) I — > (r, 9,z,u-\- s), 



gi{s) : {r,9,z,u) 



4r' 



{2s-rsm29Y 



1/2 



arctan ( -(1 + tan 
\r 

95{s) ■■ {r,9,z,u) 

(l(T 



2 6')i/2 -tan6l) ,z,i}j, 



(l+tan2 61)1/2 A Vs(l+tan2 6i)i/2-r 



r tan 9 



2\ 1/2 



arctan I 



r tan^ 



(l+tan2 0)i/2 + r 
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geis) : {r,6,z,u)i — > cos^ {s + C2) + sin^ {s + C2)j , 



{df - Ci)i/2 cos(s + C2) 
((ci - d2)cos2(c2))' tenh-i (Iii^cos(c2)) 



(d?-Ci)l/2cos(c2) 

~2" 



■((ci -d2)sin(2s + 2c2)) 
57(5) : {r,e,z,u)t — > ^icos^is + C3) + dlsin^{s + C3)j , 



(d2-ci)i/2 C0S(S + C3) 
((ci - di)cos2(c3))' tenh-i ((i%^cos(c3)) 

(di-Ci)l/2cOs(c3) 



— (^(ci -d^)sin(2s + 2c3)j 



where we assume that Ci — + , Ci = arctan(z(r^ cos^ t) ^^"^^ C3 = — s — 
arctan(2;(r^ sin^ t)^-'^/^), d\—r sint and , d2 — r cost. 

It is "well-known that for each one-parameter subgroup of the full symmetry 
group of a system there will correspond a family of solutions that are called 
invariant solutions [101 HH [E] ■ In section four, we introduce a general solution 
of CHE with respect to its functionally independent invariants. 



Theorem 3. If u — f(r, 6, z) he a solution of (QJ), then are the below functions 

gi{s)- f{r,0,z)^ f{r,e + s,z), 
g2{s)- f{r,e,z)^ f{r,e,z + s), 
gsis)- f{r,e,z)^ f{r,e,z)^s, 
5.(s)-/(r,(?,z)-/(g,(s).(r,(?,z)), 

when gi{s) ■ (r, 9, z) means the first three components of gi{s) ■ (r, 0, z, u) and s 
is an arbitrary parameter. 
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3 Invariants and General Solutions 



In this section, we study the invariant solutions of various combinations of vector 
fields ([6]) based on the following method: 

the invariants u — I(r,9, z,u) of one-parameter group with infinitesimal 
generators in the form of ^ satisfy the linear homogeneous partial differential 
equations of first order: 

v[I] = 0. 

The solutions of the latter are found by the method of characteristics (See and 
[TT] for details). So we can replace the latest Eq. by the following characteristic 
system of ordinary differential equations 

dr dO dz du 

6 " 6 " 6 " 7' 

By solving the Eqs. ^ of differential generators ([5]), we (locally) find the 
following general solutions 

(ci sinff + ca cos6')r2 +C6r 

Ii{r,e,z,u) = 2z + - ■ ' ^ ^=^1' 

(ci z + C3) sm u + (c2Z + C4) cos 

l2{r,9,z,u) = z — (ci sin 6* + C2 cos^) In f (cn BY(V— C2, Vca + fc^) 

Cl3 ^ 

+C12 BJ(V-C2, VC3 + fc^))(c7C9 exp(0^6' + y/c^ z) 
+C7C10 exp{y^ 9 - y/c^z) + C8C9 ex.p{-y/c^9 + ^/c^ z) 
+C8C10 exp{-y/(^9 - y/c^z)) + ci3 uj = d2, 
I:j{r,9,z,u) = — — arctan(a) (2 C5 ^""'^ (C2C3 — CiC4)r'^ + 2c6 r) 

A-^ ( ln(C) - In ( ^— -) ) (c^ z + C1C3 + C2C4 + c? z) 

V V COS 9 + 1/ / 

^ 2 A-^ /cos 6* — 1\ , , , 

+2r A arctan — (C2C3 — C1C4) — d^. 

\ sm 9 / 

when we assumed that 

A — c^z'^ + 2 C2C4 z + C4 + 2^ + 2 C1C3 z + C3 

D/'o 2 22i22 22o 2\i 

= (— 2 C1C3 z — C3 — z +C5r — C2 z — 2C2C4Z — €4)2 

C = (ci z COS 6* + C3 COS 9 — C4 sin 9 — C2 z sin 6* + C5 r) 

cos ^ + 1 

(ciz — ciz COS + C3 — C3 cos 9 — csr + c^r cos 61 + C2Z sin 9 + C4 sinfl). 



B sin 

and di s (i = 1, 2, 3) are some constants. 

The functions /i , /a , -^3 form a complete set of functionally independent in- 
variants of one-parameter group generated by ([8|) (see [lO]). 
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According to the theorem of section 4.3.3 of [B], the derived invariants ^ 
as independent first integrals of the characteristic system of the infinitesimal 
generators (O, provide the general solution 

S{r, 6, z, u) := fi{Ii{r, 9, z, u), hir, 9, z, u), h{r, 9, z, u)), 

with an arbitrary function fj,. 

4 Classification of subalgebras 

It is well-known that the problem of classifying invariant solutions is equivalent 
to the problem of classifying subgroups of the full symmetry group under con- 
jugation in which itself is equivalent to determining all conjugate subalgebras 
[TOl [T2] . The latter problem, tends to determine a list (that is called an optimal 
system) of conjugacy inequivalent subalgebras with the property that any other 
subalgebra is equivalent to a unique member of the list under some element of 
the adjoint representation i.e. [} Ad(g) f) for some g of a considered Lie group. 

For finding one-dimensional subalgebras, this classification problem is essen- 
tially the same as the problem of classifying the orbits of the adjoint represen- 
tation. 

The adjoint action is given by the Lie series 

Ad(exp(sXO)Xj = X, ~ s[X,, X,] + y [X„ [X„ X,]] - • • • , 

where s is a parameter and i,j ~ l,---,7. Let Afi(s) : g — > defined by 
X ^ Ad(exp(s Ari))X be the adjoint representation on g (for i = 1, • • • , 7 and 
small parameter s) and also Mi(s) be its matrix of representations with respect 
to the basis {Xi, • • • , Xy}, by an abuse of convention. After straightforwardly 
computations, we find that 



Mi(s) = 



M2{s) = 



' 1 





























1 





























1 





























coss 




sin s 





















— sin s 


cos s 































cos s 


— sin s 





















sin s 
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" 1 





























1 





























1 





























1 





























1 


























-s 


1 




















—s 








1 








Id7, 





















8 



M5{s) = 



1 








—s 







" 
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COS s 
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— sin s 























cos s 
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In the matrices Mi,Mq and My, there are minor matrices which indicate 3- 
dimensional rotations, that are related to the section so(3) of the Lie structure 



Theorem 4. An optimal system of one- dimensional Lie subalgebras of is 
provided by those generated by 

1) Xi, 10) aXi + hXQ, 

2) X2, 11) aX^ + hX-j, 

3) Xq, 12) X:i + aXi+hX2, 

4) X-!, 13) X-i + aXi + hXfi, 

5) X3 + aXi, U) aXx + bX2 + cX^, 

6) X^i + aXi, 15) aX2 + bX5 + cXT, 

7) Xs + aXe, 16) aXi + bX^ + cX^, 

8) Xa + aXy, 17) X:i + a Xi + b X2 + cX^, 

9) aXi + bX2, 
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where the coefficients a, b and c are arbitrary constants. 

Proof: Let X = X^J^;^ a, Xj be an element of g, then we see that 

M-j{s7) o ■ ■ ■ o M2{s2) o Mi{si) : X i — > (^ai cos Sq cos s-j -\ + a7sins6)-^i 

+ ( 02 COS sq cos sr + tti COS sq sin 57 + 05 sin sq ) X2 



+ \a\ (sin si sin S7 — cos s\ sin se cos S7) H \-aj cos s\ cos s%j X7. 

For classifying one-dimensional Lie subalgcbras of CHE, we plan the following 
cases that in each case, by acting a finite number of the adjoint representations 
Mi [i = l,---,7) on X, by proper selection of parameters Si in each stage, 
we gradually try to make the coefficients of X vanish and to acquire the most 
simple form of X. 

• At first, wo suppose that 03 7^ and scaling if necessary, we can assume 
that as = 1. This assumption suggests Case 1,- ■ ■, Case 10. 

Case 1. If 05 7^ 0, then we act on X by Ad(cxp(— arctan(a4/a5) Xi)) and 
hence we can make the coefficient of X4 vanish. Then we tend to the new form 

X' = aiXi+ 02 X2 + X3 + a's X5 + 03X6 + a'^ X-j 

for certain scalars 05 and Og, 07 depending on 04, 05 and 04, 05, ag, 07 resp. By 
acting Ad(exp(^ -^4)) on X' when we suppose that a2 ^ 0, we tend to 

X" = a'l Xi + 02 X2 + X3 + q/, A5 + o!^ Xe 

where the coefficient of X5 is vanished for scalar a[ depending to ai, a2, Og, a^. 
Also be the action of Ad(exp(— (ag/a2) X^)) on X" the coefficient of Xq will be 
zero: 

X'" = a J Xi + a2 X2 + X3 + Og Xr,. 

At this stage, by acting adjoint representations Mi{ai) on X'" , we find that no 
more simplification of X'" is possible. Thus each of a'^ and ag are arbitrary. 

Therefore, any one-dimensional Lie subalgebra generated by X with 02, ag 7^ 
(by knowing that 03 7^ 0) is equivalent to the Lie subalgebra spanned by 

a'l Xi + 02 X2 + Xs + Xg, 

which is equal to section 1 7) of the theorem. 

Each of the following cases is prepared by a similar method to Case 1 and 
by eliminating unnecessary details, we just give the conditions. 

Case 2. If as 7^ and 02 = 0, then we can make the coefficients of X6,Xi 
vanish resp. by M2, M4 and then X is reduced to 11). 



10 



Case 3. In the case which = along with a2,a4 ^ 0, we act on X by 
M4, M5, M5, M7 to cancel the coefficients of Xt,Xq, Xi,X2 resp., then the sim- 
plest possible form of X is equal to the part 6) of the theorem. 

Case 4- When we change the condition 04 ^ in Case 3. to 04 = 0, by applying 
AI4, Ms the coefEcients of X^, Xq vanish and no more simplification is possible. 

Therefore X is reduced to 12). 

Case 5. When a2 = 05 = and 04,07 ^ 0, then by applying M4,Mq,M2 one 
can make the coefficients of X5, Xi, X-j vanish and hence X turns to 13). 

Case 6. If in Case 5 we only change 04 ^ to 04 = along with ae 7^ 0, then 

we can act on X by resp. M4, Afg, ^1 to cancel the coefficients X5, X\, X7. 
Hence we find that form of X is in the form of 7) . 

Case 7. In the latter case, if we suppose that oe = 0, by acting M4, Mq the 
coefficients of resp. X5, X\ are zero and X is equal to the part 8) of the theorem. 

Case 8. Let a2 = 05 = 07 = and 04 ^ 0. we can make the coefficients of 
X5, Xi resp. by M4, M5 vanish and then reduce X to the section 13). 

Case 9. When a2 = 04 = = a-j- = and 7^ 0, by acting resp. M4, M7 on 
X, one can make the coefficients of X^, Xi vanish. With these conditions, X 
gives the form of 7). 

Case 10. When a2 = 04 = 05 = ae = 07 = 0, by the action of M4 on X and 
vanishing X5, X is reduced to 5). 

• At the second, we suppose that as = and consider Case 11,- ■ ■, Case 20 

prepared by this assumption. 

Case 11. Let 04,05 ^ 0. We can act by M2, M7, M4 to cancel the coefficients 
of X7, X2, Xi resp. then we will have the reduced form 16) of X. 

Case 12. By changing 05 ^ of the latter case to be zero and applying M2, 
M7, M5 on X, the coefficients of X^, X2, Xi resp. will vanish and then X is 
reduced to 10). 

Case 13. Let 04 = 0, 07 ^ and 05 ^ (by following the process and finding 
the exact generating coefficients, in fact, we assume that a'^ = 0102 + 0507 be 
non-zero). The coefficients of Xi, Xe by acting Ma, M2 will be zero and no 
further simplification is possible. Hence X is reduced to 15). 

Case 14. Suppose that 04 = Og = 0, 07 ^ 0, and 02 ^ (in fact, in the method 
of construction, 02 = —0207 -|- 0105 ^ 0). Applying MQ,Mi,Ms resp. on X 
make the coefficients of Xi, X7, Xq vanish. In this case X is reduced to 2). 

Case 15. In the last case, we let a'2 be zero and apply Mg, Mi, Mi resp. to 
cancel the coefficients of Xi, X2, X^. Then X will be in the form of .^). 

Case 16. If O2 = 04 = Og = 67 = then by acting Mg, Mi resp. to cancel the 
coefficients of Xi, Xq, we lead to the simplest case 3) of X. 
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Case 17. Let 04 = ay = and ^ 0. By acting M2 on X then the coefficient 
of Xg will be zero and hence X is reduced to 14)- 

Case 18. If a4 = 05 = 07 = and 02 ^ 0, then by acting A/5 to cancel the 
coefficient of Xg, we lead to the simplest case 9) oi X. 

Case 19. When 02 = 04 — 05 = 07 = and ag 7^ 0, then we can make the 
coefficients Xi resp. by Mj vanish. Then X is reduced to 3). 

Case 18. Finally, if in the last case ag — 0, no further simplification is possible 
and then X is reduced to i ) . 

There is not any more possible case for studying and the proof is complete. A 

For finding optimal systems of higher dimensional Lie subalgebras (equiva- 
lently Lie subgroups) , we refer the reader to [12] or [7] . 
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